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Abstract 

An analysis of radiating perfect fluid models with asymptotically AdS bound- 
ary conditions is presented. Such scenarios consist of a spherical gas of radiation 
(a "star" ) localised near the centre of the spacetime due to the confining nature 
of the AdS potential. We consider the variation of the total mass of the star as 
a function of the central density, and observe that for large enough dimension- 
ality, the mass increases monotonically with the density. However in the lower 
dimensional cases, oscillations appear, indicating that the perfect fluid model 
of the star is becoming unrealistic. We find the critical dimension separating 
these two regimes to be eleven. 

1 Introduction 

There are numerous effects in physics which are dependent on the dimension of 
the spacetime in vi^hich they live, and concepts which have mathematically simple 
models when restricted to one or two dimensions may become intractably complex 
as the dimensionality increases. Furthermore, in some situations there is a sharp 
contrast in the behaviour of the system when the number of dimensions passes 
above or below a certain "critical dimension" . In the field of general relativity there 
are already some instances of such dimension dependent phenomena; for example, 
the Gregory-Laflamme instability [H 12 of black strings, and the work of Belinsky, 
Khalatnikov and Lifshitz (BKL) [H] and its extensions [HOE], where the dynamics 
of a spacetime in the vicinity of a cosmological singularity were studied. In the latter 
case they found that the general behaviour of the relevant Einstein solutions changed 
from "chaotic" in the low dimensional cases to non-chaotic in higher dimensions. It 
is this extra complexity and appearance of a critical dimension which is discovered 
in the work presented here. 
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The issue of dimensionality has become of even greater importance over the past 
few decades, with the theories of hidden dimensions first postulated by Kaluza [7] 
resurfacing in the quest for unification and a Theory of Everything. Originally it 
was hoped one extra dimension would be sufficient; with the advent of string theory 
and its subsequent developments this was then expanded to twenty six in the late 
1960's as a consistency requirement for bosonic strings, before being reduced back 
down to ten with the introduction of supersymmetry in the 1980's. The idea of 
holography also has its roots in dimensionality, and its most famous form, that of 
Maldacena [Hj, states a relation between string theory on a five-dimensional anti-de 
Sitter space, and a four-dimensional conformal field theory. Much work (see [9] for a 
review) has followed investigating spacetimes and conformal field theories of varying 
dimension and complexity. 

Here we present an analysis into the stability of radiating perfect fluid spheres 
in an asymptotically anti-de Sitter spacetime. Due to the confining nature of the 
AdS potential, the spheres of radiation are self-gravitating, and thus we shall refer 
to them as "stars" in much of the following, although this is mainly used as a 
conveniently brief label, as it is only a toy-model approximation to a star at best. 

We begin by giving the equations for a such a model in d-dimensions, before 
analysing the behaviour of the star's total mass. We consider the variation of the 
total mass as a function of the central density, and observe that for large enough 
dimensionality, the mass increases monotonically with the density. However in the 
lower dimensional cases, oscillations appear (this was originally noted in the d = 5 
case in |10]). indicating that the perfect fluid model of the star is becoming unre- 
alistic. We numerically find the critical dimension separating these two regimes to 
be 11.0 (to three significant figures), and give an explicit relation, ()lip . between the 
spacetime dimension d and the "saturation density" pc, see section 13.11 The exis- 
tence of local maxima at critical central densities (saturation points) in the lower 
dimensional cases indicates the appearance of instabilities in the model of the star, 
and point to unrealistic effects developing for pQ > pc- We also provide a numeri- 
cal analysis of the behaviour at large central density, in particular the self-similar 
behaviour that appears in dimension d < dc, several parameters of our numerical 
model are then also determined analytically from a dynamical systems analysis of 
the behaviour, where we consider the expansion about a fixed point of the zero- 
cosmological constant solutionis 

The outline of the paper is as follows: we begin with a brief recap on perfect 
fluid models in general dimension in section [2l In section [3] we introduce the total 
mass as a function of the central density, analyse the progression of the saturation 
point with increasing dimension, and present the best fit formula, (jlip . which gives 
a value for the critical dimension. We also present further numerical results for 

^This dynamical systems analysis (section Q was suggested after correspondence with 
V. Vaganov, who also considered the behaviour of self-gravitating radiation in AdSd in [11]. Work 
simultaneously conducted by P. H. Chavanis also found the critical dimension described here, via 
an alternative route, in his comprehensive study of relativistic stars with a linear equation of state, 
[12j (see the note following the discussion for more details). 
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the behaviour of the total mass at large central densities. In section H] we give a 
dynamical systems analysis, following the methods in [11] and [13], which yields 
analytical expressions for several of the numerical results of the previous section. 
Finally, we conclude with a discussion of the results in section [5j 



2 Perfect fluid models 

Consider a general static, spherically symmetric d-dimensional AdS spacetime with 
metric: 

ds^ = -k{r)df + h{r)dr'^ + r'^d^l_^ (1) 

By considering a perfect fluid of a gas of radiation, one can obtain implicit expres- 
sions for k{r) and h[r) for a simple model of a "star" geometry. For a perfect fluid 
we have that the stress tensor is of the form: 

Tab = p{r)UaUb + P{r){gab + UaUb) (2) 

where u"" is the d- velocity of a co-moving gas, and upon which we impose the further 
constraint that the matter be purely radiating; this sets p{r) = (d — l)P{r) as it 
requires that Tab be traceless. One obtains the required metric by solving Einstein's 
equation: Gab + ^dab = ^irGciTab, with the above stress tensor and a negative 
cosmological constant, as follows. The relevant components of Einstein's equations 
in general dimension d are given by: 



fd-2\ {d-3){k{r)-k{r)h{r))+rk'{r) f p{r) , (d - l){d - 2) \ ^ 

(3) 

(4) 

where we have used that A = -{d - l){d - 2)/{2R^), P{r) = p{r)/{d - 1), and set 
SirGd = 1 for convenienc^. We can infer the form of h{r) from ([3]), as the k{r) 
dependence cancels, and we find that h{r) is given by: 



where m(r) is a mass functiorH related to the density 



via: 



^In the numerical results presented shortly we also set _R = 1; we include it here to ease com- 
parison with the dynamical systems analysis of the A = {R = oo) case given in section O 

^Note that the mass function used here is a rescaling of the actual mass; a constant factor 
from the integral over the angular directions does not appear in our definition of m(r) due to our 
definition of h{r). 
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In order to specify a form for k{r), we recall the energy-momentum conserva- 
tion equation, V^T^'^ = 0, which for a general perfect fluid without the radiation 
condition gives: 

P'{r) + ^^{P{r)+pir))=0 (7) 
which can be re-arranged to give 

in the radiation case, where we have introduced poo, which is the leading coefficient 
of p{r) at large r, and is given by poo ~ p{r)r'^ as r — > oo. Substituting h{r) from 
^ into dS]) and eliminating k'{r)/k{r) using ([7]) then gives an equation in terms of 
m{r), p{r) and p'{r),. 



(d-S) / 1 \ 2p'{r) 2p{r) + {d-lf{d-2)/R^ 



(9) 



V l + 'P^'^d (d-l)(d-2)(l + f 

which couples with our equation for m'{r): 

m'{r) = p{r) r'^''^ (10) 

to give a pair of ODEs. For specified dimension d, these allow the geometry of the 
spacetime to be numerically generated when they are combined with the relevant 
boundary conditions: m(0) = and p(0) = po- The condition p(0) = po specifies 
the central density of the gas, and we have that for fixed R, po is the single free 
parameter of the system (pure AdS is recovered when pQ = 0). 



3 Total mass as a function of central density 

Whilst mathematically one can work with this perfect fluid setup in any number 
of dimensions (including non- integer ones), one would also like to consider the ap- 
propriateness of doing so, given that we wish to use the geometry as the setup for 
a toy model of a star. In other words, is there any significant change in behaviour 
as the dimensionality of the model is altered. A particular quantity of interest in 
analysing the stability of the model is the total mass M of the star, and as we have 
just seen, the mass and density profiles of our gas of radiation are determined by a 
single parameter: the central density of the gas, pQ. 

To avoid possible instabilities such as those considered in the asymptotically 
flat case (in four dimensions) in [T3], one would expect the total mass to increase 
monotonically with pQ. One could also expect the total mass to be bounded from 
above by some maximum value, analogous to the Ad asymptotically flat case where 
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Fig. 1: Total mass vs density for the radiating perfect fluid model in various dimensions, from 



d = 4 (top curve) through to d = 12 (bottom curve). The saturation point for each 
dimension is indicated by the red dots; these correspond to the maximum value of the 
total mass in the relevant dimension, at the critical density pc (see section [3?l]l . For d 
large, there is no local maximum and hence no finite saturation point; in these cases, 
the maximum total mass is given by the asymptotic value, rjd- 



for a fixed size Rstar, the maximum possible mass sucli a star can have is given by 
Mmax = 4:Rstar/9, a result found by Buchdahl in 1959 [15] Fl 

In our scenarios the total mass is indeed bounded from above, however, this 
maximum is not always the asymptotic value of the total mass at large density. 
Although we observe that as po ~^ oo we have M{pq) r/^, where r]d is some finite 
constant dependent on the dimension d (see section [32]below for more details), what 
we do not find in all cases is the total mass approaching this constant monotonically, 
see figure [TJ When the dimensionality is low, there are sizable oscillations about the 
final value rfd before the curve settles down (see figure [2]), as was noted in the d = 5 
case in [TO] , and in other similar scenarios, e.g. [IBj, and the star's maximum mass 
is given by some value greater than r/^. As the dimension is increased, however, 
these oscillations become less pronounced, and for d sufficiently high they disappear 
altogether, see figure El 

Ideally one would like to analytically determine the dependence of the shape 
of the curve on both the dimension d and the density pQ, however, due to the 
complexity of the equations, the exact behaviour must be computed numerically. 
One can nevertheless use this data to construct models of the various features of 
the star's behaviour: for example, in section 13.21 below, we give an analysis of the 

^For further edification, note that this equaUty coincides with the Buchdahl-Bondi limit, usually 
written in the form R/M = 9/4, which is the lowest radius to which the Schwarzschild geometry 
can be embedded (in Euclidean space). For more detail see e.g. [TBI [TT] . 
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Fig. 2: The oscillations in the total mass M: as the central density po is increased, M does 
not simply increase monotonically towards its final value rjd. Instead, it reaches a larger 
maximum before undergoing damped oscillations towards Note that the amplitude 
of the oscillations becomes smaller as the dimensionality d is increased. 
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Fig. 3: For larger d, there are no oscillations in the total mass: M is now a monotonic function 
of the central density po. and its maximum is also its asymptotic value as po — * oo, 
namely r]d- 
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total mass at large (where it approaches a constant, dependent on d) in different 
dimensions. 

One particularly interesting feature is the appearance of the turning points in 
the total mass seen in figure [21 and specifically the locations of the local maxima 
in different dimensions. One can see from the figures that as the dimensionality is 
increased, the appearance of the first maximum moves to larger pq; by analysing this 
progression one can obtain a remarkably simple relation which immediately gives a 
value for the critical dimension, above which the oscillations do not exist, and hence 
the total mass is a monotonic function of po- 

3.1 A critical dimension 

The saturation point, po which we define as being the location of the first local 
maximum when increasing pQ, can be seen to progress to larger and larger po as the 
dimension d is increased, see figure [H What we wish to determine is whether this 
saturation point appears for all dimension d (for sufficiently large po), or whether 
there is a cut-off dimension, dc, such that for larger d, there is no local maximum 
and hence no saturation point. Figure H] shows how the saturation point varies with 
dimension; numerical analysis then finds (to 3 significant figures) that this behaviour 
is given by the following model: 



which gives a critical dimension dc = 11. What is perhaps rather surprising is 
the simplicity of pi|) : not only do we have a critical dimension appearing so clearly, 
the overall dependence on d is remarkably simple, and the co-efficient of the linear 
term appears to be exactly one half. 

As we shall see in section [U the value of the critical dimension can also be 
determined by an analytical consideration of the radiating perfect fluid system with 
zero cosmological constant (i.e. in the limit R oo). Although such a solution is 
singular at r = 0, and has infinite mass, by confining the radiation to finite sized box 
one can obtain finite mass solutions. The features determined in this configuration 
can be related to equivalent behaviour in the asymptotically anti-de Sitter case 
(where the (finite) mass is confined by the AdS potential), and indeed exact values 
for certain parameters can also be computed. This same analysis is not restricted to 
the star geometries considered here, it can be used with any linear equation of state 
|llj . or even more generally [13]. Before giving the analysis for our case of perfect 
fluid radiation, however, we firstly present further numerical results. 

^As mentioned earlier, correspondence with V. Vaganov and P. H. Chavanis suggested that the 
critical dimension in the radiating perfect fluid case is very close to (but not exactly) eleven, and 
this is indeed the case as we see in the dynamical systems analysis approach in section |4l where 
we obtain a value complementary to the numerical estimate of 11.0 given here. Interestingly, the 
exact value of dc = 11 appears in the case of Newtonian isothermal spheres, as noticed by Sire and 
Chavanis in 2002 [19]. 
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Fig. 4: The change in the saturation point pc with increasing dimension d. The points plotted 
are the calculated values for the saturation point for the star model in the corresponding 
dimension, the red best fit line is the curve given by (|lip . The divergent behaviour as 
d approaches eleven indicates that for d > 11 there is no saturation point, and hence 
no apparent instability in the perfect fluid model of the "star". 
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Fig. 5: The plot on the left shows rjd for various dimensions, with the approximation given in 



shown in red. The data points are all at integer values for the dimension, with 
the addition of points at d = 3.1, 3.2, 3.5 to highlight the behaviour of the curve 
at low d. The righthand plot shows the behaviour of Af™"^; this is identical to that 
of 77(j for d > 11, however for < 11, the maximum is given by the value of the total 
mass at the saturation point, pc- The best fit approximation (red curve) for each is 
simple in terms of its d dependence, and provides a good fit over a large range of d. 



3.2 Total mass at large po 

In addition to considering the variation of the saturation point for the star with 
dimension, one can also investigate the asymptotic behaviour of M as pQ becomes 
large. As mentioned in section [3l at large po, the value of the total mass tends to 
a constant, rjd, which is then only dependent on the dimension; the value of this 
constant decreases as d increases. The values are plotted in figure [5] and despite 
the complicated appearance of the plot, a remarkably close fit for all dimensions is 
given by: 



which is also shown in the figure. Checks show that the function continues to 
give accurate predictions for larger d, and although there is perhaps slightly more 
complicated behaviour for d ~ 4, we do not attempt to investigate this further her^; 
despite the relative compactness of the expression, there is little intuitive origin for 
any of the constants involved. Nonetheless, it is impressive that the behaviour of 
the mass at large po can be expressed in such simple powers of the dimension. 

One can perform a similar analysis of the behaviour of the maximum value of 
the total mass as the dimension increases; the results are also shown in figure [5j For 
d > 11, the maximum total mass corresponds to the asymptotic value, ijd, however, 
for lower dimension, the maximum is given by the mass at the saturation point. A 

^Whilst in the dynamical systems analysis (section the d = 3 case (where we have r/a — 1) 
needs considering separately, as there is different asymptotic behaviour involved due to the non- 
dynamical nature of gravity in such a scenario (see 20 for example) , we find that we can include it 
both here (in the analysis of rjd) and also in our earlier result for the critical dimension (see section 



% « 0.716 + exp [9.85 - 3.72 d] - exp -0.603 
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0.371 


8.44 
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1.03 
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Table 1: Numerical estimates for ad, /3d, IJ-d and Vd (to three significant figures) for the model 
of the total mass given in ([T4|) . 
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good fit to the curve is given by: 

20.5 



^max _ Q_7J2 + exp [2.74 - 1.07 d] - exp 



-0.592 



d 



which differs (significantly) from (I12p only in the second term; this was to be ex- 
pected, as the curves differ only at low d. Whilst their is little apparent significance 
about the values of the numerical constants involved in the expression, we have 
again produced a fit with a relatively simple dependence on d which gives accurate 
predictions for M™"^ over a large range of dimensions. 

Although the form of the fit used in equations (fT2|) and (fT3l) was chosen primarily 
because it gave such a close fit to the data, it would be interesting to study possible 
reasons for expecting the observed dependence on dimension, although we do not 
pursue this further here. What we do now examine in more detail, is the oscillatory 
behaviour, which can be considered both numerically and analytically. 

3.3 Self-similarity analysis for d < 11 

Another interesting feature of the plots of the total mass seen in figures [1] and [2] is 
the self-similarity exhibited by the oscillatory behaviour as pQ ^ oo. A numerical 
analysis of the periodicity and damping of the oscillations seen for 3 < d < 11 leads 
us to propose the following model for the total mass: 

Md{po) ^rjd + Od exp [-(3d log(/9o)] cos [pd - Vd log(po)] (14) 

which gives a good approximation for the behaviour in the region po > Pc- In (I14p . 
rjd is the asymptotic value of the mass discussed above, and the four parameters 
Oid, Pd, Pd and Vd are constants dependent only on the dimension d. Approximate 
values for these constants for d = 3.1, 4, 5, 6 and 7 are given in Table [TJ 

Although the values given in table [T] are only approximate, we nonetheless see 
interesting dependencies on d emerging. For example, Pd appears to increase roughly 
linearly with dimension (Pd ~ 0.22d — 0.51), as do ad and pd for d > 4. We will 
see below in the dynamical systems analysis how this linear behaviour of Pd on d is 
only an approximation to the true behaviour, and the same analysis also provides 
exact values for the parameter i^d- This analytical analysis also confirms the form 
of the fit used in the numerical approximation derived above. 
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4 Dynamical systems analysis 

By considering the behaviour of the system of coupled ODEs given in section [2] 
in the Umit i? — > oo, we can obtain analytical results for some of the interesting 
features of the radiating perfect fluid star geometries described above. The analysis 
presented here follows that detailed in both p[T] and [13], where it is given in more 
general settings; by focusing on the radiation case (where p{r) = (d — l)P{r)) we can 
give a good explanation of why the numerical behaviour seen above is so, without 
excessive over-complication. 

The basic idea is to rewrite the equations for p'{r) and m'{r) in terms of di- 
mensionless (compact) variables and perform an analysis of the fixed points. The 
corresponding eigenvalues and eigenvectors obtained by linearising about these fixed 
points give a complete description of the nearby behaviour (Hartman-Grobman the- 
orem, |21j ) on the new state space, which can then be translated back to the physical 
picture by inverting the transformations given below. Interestingly, for the perfect 
fluid stars, the dependence of the total mass (as well as other quantities, e.g. the 
entropy) on the central density, po, is governed by the behaviour around (and hence 
the eigenvalues of) a single fixed point. Specifically for our case we will see how 
this gives both an exact value for the critical dimension dc, and a clear analytical 
explanation for the observed behaviour in the two regimes d < dc (oscillatory) and 
d > dc (monotonically increasing). We will also obtain expressions for the (3^ and 
I'd parameters introduced earlier. 

To proceed, we thus set i? = oo, and our equations Q and (fTOj) become: 



where SirGa has again been set equal to one. Note that we do not include the d = 3 
scenario here as it is a special case (due to the non-dynamical nature of gravity). 
We can now introduce the dimensionless variables: 




)) 



(15) 




(16) 



u = 




(17) 



and 




(18) 



which allow equations ()15p and (I16p to be rewritten in the form: 




(19) 




(20) 
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where we have also introduced the new independent variable ^ = ln(r). For the case 
of positive mass and density we're considering here, both u and v are greater than 
zero (for non-zero r), and we make a final change of variables to the bounded U and 
V defined hM 

U = -^, V = ^ (21) 



1 + t; 



l + u 

which gives the system of equations: 

(22) 
(23) 



dV 



V{1 -V){3-d+{d- 2)U) ( 1 + , , , 

\ \d — 1 



where we have also introduced the independent variable A, defined by: 



d\ 



{l-U){l-V) 



{l + u){l + v) 



(24) 



The fixed points of the system ( (j22p and (j23p ) are calculated in the usual fashion, 
by setting both dU/d\ and dV/dX to zero and solving for U and V\ there are six in 
total, with eigenvalues then obtained from 



d_ 
dX 



U 
V 



a (du\ 

dU \d\J 
d (dV\ 
aU \d\J 



a (du\ 
ay \d\j 

A. (dV\ 
dV \d\J 



fp 



U-Ufp 
V-VfP 



(25) 



where the matrix components are evaluated at the particular fixed point under 
consideration. A table of such eigenvalues is given in pT], where they are labelled 
Ti, . . . ,Tq; we do not list them all again here, however, as orbits in the interior 
of the state space [0, 1]^ originate from either T2 or T4 and converge to the fixed 
point T3 (as is shown in [E]). This fixed point corresponds to the singular self- 
similar solution given by equation (2.14) of [llj . and due to the scale invariance 
of the system one can consider the entire set of (positive mass) solutions as being 
represented by a single orbit from T2 to T3 (this is true for any linear equation of 
state, P{r) = qp{r)). 

Although one cannot write an analytic expression for this orbit, one can obtain 
approximations by linearising about the fixed points. As discussed briefly earlier, as 
this zero-cosmological constant solution is singular, in order to produce finite mass 
solutions the radiation must be confined to an (unphysical) box; the two fixed points 
T2 and T3 thus represent solutions with po = and in the limit /?o — >■ 00 respectively. 
The behaviour described by the linearisation about T3 then reveals aspects of the 
large po limit of the radiating stars (where the confining AdS potential results in 



^Although the range of both U and V is defined as being (0, 1), in order to perform the fixed point 
analysis of the asymptotic behaviour, it is necessary that the boundary points also be included; 
this requires the system given by (|22|l and (|23p be on [0, f]^, which is manifestly so. 
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the finite mass solutions without the need for any unphysical box), exactly what we 
analysed numerically in sections 13.21 and 13.31 This linearisation gives an explanation 
for the existence of a critical dimension and the differing behaviour seen in higher 
and lower dimensions, including quantitative expressions for dc and the and i^d 
parameters of Table [H as we shall now show. 

Fixed point T3 corresponds to the following values of U and V: 



TT T/ 2{d-lf 

^^^ = d^' ^^^ = 2-4d + (d-l)d2 (26) 



and has eigenvalues: 



_ d(d-3) ( ^ ^ , / (d- 12)^2 + 13 0; -18 
^3 - 2-Ad+(d-l)d'^ V d-2 




(27) 

where we denote the coefficient and observe that it is strictly positive for d> 2>. 
These eigenvalues govern the behaviour of the solution, and we immediately see that 
there are two distinct regimes; one where the expression inside the square root is 
negative, corresponding to the oscillatory behaviour seen in figure [21 and one where 
the expression is positive, resulting in the monotonic behaviour seen in figure [sj^ 
We thus obtain a value for the critical dimension given by the solution to: 

{dc - 12)dl + 13 4 - 18 = (28) 

which yields dc = 10.964... , complementary to the value of dc = 11.0 obtained 
numerically, although with the significance of being non-integer rather than exactly 
11; interestingly for any linear equation of state the value of dc is always in the range 
10 < 4 < 11, see dl]. 

We can relate the asymptotic behaviour obtained from the state space picture to 
the physical quantities of mass and density via several auxiliary equations to those 
given above, specifically: 

dv nun 

_ = ^ = u{l-V)m (29) 

and 

dp Vd f U \ 

Given expressions for U and V in terms of A (as obtained from an analysis of the 
behaviour around the fixed points, say), one can integrate the above to determine 



The fact that kd > ensures that the fixed point Ta is a stable focus for the oscillatory behaviour 

3 



in the d < dc case; for d > dc, we have that is strictly less than zero, and hence acts as a stable 



node. 
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corresponding expressions for the mass, radius and density in terms of A. There is, 
however, a simple way to see the dependence of the total mass on the central 
density po, which also reveals the origin of the (3d and i^d parameters of our numerical 
model in section 13. 3[ 

Focusing then on the case where d < dc, how does the imaginary term in ()27p 
lead to the (self-similar) oscillatory behaviour manifest in the total mass at large po? 
This can be seen directly from the linearisation about T3, where we observe similar 
oscillations in the expressions for U{\) and V{\) (see below); as mentioned above, 
it is the behaviour around T3 that governs the behaviour of the physical quantities 
in the large po limit. By considering the behaviour of U{X) and V{X) in terms of 
Pq, we can extract the coefficients which should then match those in (jl4p (as argued 
more fully in [E]). The solutions of ([22]) and ()23p in the large A limit (i.e. about 
the fixed point Ts) can be expressed as: 



Re{UiX)) = C/T3 +exp(-(ci- l)fc,A)cos ( kA\l ^^ ^^^^ 



Re{V{X)) = Vt, + exp {-{d - l)kA) cos (^kA^ l^_{±_^^}d^^^ ^32) 

where we have only kept the real angular term as we are only interested in the 
period of the oscillations (i^d) a-nd the coefficient of the damping (/3rf); the extra 
factors (namely ad and fid) cannot be extracted directly from this analysis If] 

For sufficiently high density stars (i.e. with large po), we have A cx | j^^^ log(/3o), 
and we thus obtain: 

_ {d - l)kd _ d 1 _ 3 

2{l-UT,)il-VT,)~ 4^ 2d 4 ^"^"^^ 

and 



kd /I8- (d- 12)^2 _ 13 d 



2{1-Ut,){1-Vt,)\ d-2 
= ^V((i-2)(18- (d- 12)d2 _ 13 d) (34) 

which give the values shown in Tabled provided as a comparison to the numerical 
estimates obtained in section 13.31 We see that they match very closely, with any 
discrepancies most likely due to a combination of numerical imprecision in the orig- 
inal data for the mass at large po and the use of oscillations at insufficiently large 
Po for the asymptotic dependence to be totally dominant. 



^Technically, for a solution of this form one should first make a linear change of coordinates such 
that the matrix on the RHS of (|25p is in diagonal form. This only manifests itself, however, as 
extra multiplicative constants which do not affect the decay term Pd or oscillation period z/^, and 
hence can be ignored. 
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d 


Ha 




3.1 


231/1240 PS 0.186 


V695519/1240 « 0.672 


4 


3/8 = 0.375 


^/47/8 « 0.857 


5 


3/5 = 0.6 


2^/6/5 « 0.980 


6 


5/6 « 0.833 


V^13/3/2 w 1.04 


7 


15/14 w 1.07 


2^/215/5 w 1.05 



Table 2: Exact values (alongside decimal equivalents) obtained from the dynamical systems 
analysis for (3d and I'd for the model of the total mass (HH). 

5 Discussion 

What we have seen in the above analysis is firstly the appearance of a critical 
dimension (dc = 11.0) from the simple requirement of a stability condition on our 
perfect fluid model for a gas of radiation, and a numerical study of the behaviour 
in various dimensions. Whilst the appearance of oscillations in the variation of 
the total mass with the central density po in the lower dimensional cases had been 
noted before, we saw here that such oscillations do not appear to persist in the higher 
dimensional cases (see figures [2] and [3]) . Not only do the oscillations die down as the 
dimensionality is increased, by analysing the progression of the saturation point, pc, 
we find that they disappear completely for d above a certain value. Although this 
value was calculated to be 11.0 in the original numerical analysis, the dynamical 
systems approach which followed in section H] not only gave a more precise, non- 
integer value {dc = 10.964 . . . , as first derived in [H]), but also explained analytically 
why one sees a change from oscillatory to monotonic behaviour as one increases the 
dimension past this point. The section concluded by demonstrating how the f3d 
and I'd parameters of the numerical model could also be derived via this analytical 
approach. 

What is remarkable is that from a seemingly basic condition (that of monotonic- 
ity in the variation of the total mass), one arrives at such a simple relation for the 
dependence of pc on d, namely equation (fTTj) . Such simplicity could not have been 
expected given the complex nature of the initial setup, which allows the spacetimes 
in question to be generated only numerically from the coupled ODEs. The exten- 
sions [T, 'F, 6j to the BKL work on modelling a gravitational field close to a spacelike 
singularity also reveal a critical dimension of eleven In their work they found that 
the general behaviour of the relevant Einstein solutions changed from "chaotic" in 
the low dimensional cases {d < 11) to non-chaotic in higher dimensions {d > 11), in 

^"Briefly, their analysis of the setup was performed using the mixmaster model, where the dy- 
namical behaviour is governed by Kasner exponents and conditions upon them, and in which the 
evolution continues until the system reaches a stability region where the Kasner exponents remain 
constant. They observed that such a stability region could only exist for d > 11, and thus the evo- 
lution continues indefinitely for any lower number of dimensions. This has interesting consequences 
not elaborated on here, which are discussed in detail in the papers cited above. 
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much the same manner as we observe the transition from osciUatory to monotonia 
total mass behaviour in the radiating star case considered here. It is interesting 
that their work also reveals a critical dimension of eleven, and a more detailed com- 
parison of the two different scenarios (including an analysis to determine the exact 
(possibly non-integer) value of dc for their transition) may yield further insight. 

It would be interesting to see if such a result appears in other investigations into 
scenarios similar to the radiating perfect fluid model considered here. For example, 
one could examine other physical equations of state to see if they exhibit the same 
behaviour, and indeed the work of |TT] and |T2] has pursued this idea further (see 
the note below). Finally, one could also look to explain the linear d/2 term which 
appears in (jlip , and whether there is any physical explanation for why the coefficient 
should take on the value of a half. 

Note 

Unknown to the author, this phenomenon has also been simultaneously investigated 
in two other works. In [llj, Vladislav Vaganov analyses the behaviour of radiating 
perfect fluid models in d-dimensional AdS spacetimes; he notes (as we do here) that 
there is a significant change in the behaviour of the total mass for d > 11 (where it 
becomes a monotonic rather than oscillatory function of the central density), and 
demonstrates this not only for the mass but also the temperature and entropy. 

He also presents a dynamical systems analysis (based on that given in [13]) of 
the behaviour for a general linear equation of state, -P(r) = qp{r), which includes 
the radiation case. This analysis complements the numerical results presented here, 
providing an analytic derivation of the critical density, which is determined to be 
dc = 10.964..., consistent with our relation (llip . The specific analysis for the 
radiation case is given in section HI where we give not only the derivation of the 
critical density, but also demonstrate how the dynamical systems technique gives 
analytical expressions for other parameters introduced in our numerical investigation 
into the self-similar behaviour for d < dc- 

The second related paper, [12] by Pierre-Henri Chavanis, presents an in-depth 
study of the behaviour of general stars ( "isothermal spheres" ) with a linear equation 
of state in an asymptotically flat background. His results are again complementary, 
finding that there is monotonic behaviour for d > 11, in contrast to the oscillatory 
behaviour observed in lower dimensions. By asymptotic analysis he also finds the 
value for the critical dimension in the radiation case to be very close to eleven, and 
although there initially appeared to be a discrepancy between the two alternative 
calculations of the critical dimension in [11] and |12j . the latter was subsequently 
corrected to agree with the value of dc = 10.964. . . found in [TTj. His paper also 
includes a comprehensive investigation into the stability of the different regimes, 
looking at a number of alternative stellar configurations and considering the be- 
haviour of other thermodynamic parameters (entropy, temperature,. . . ), in addition 
to the mass. 
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